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The evolution of the gluon plasma produced with satura- 
tion initial conditions is calculated via transport theory for 
nuclear collisions with 0.1 < y/s < 10 A TeV. The effective 
longitudinal pressure is found to remain significantly below 
the lattice QCD pressure with these initial conditions until 
the plasma cools to near the confinement scale. The abso- 
lute value of the transverse energy per unit of rapidity and its 
dependence on beam energy is shown to provide a sensitive 
test of gluon saturation models since the fractional transverse 
energy loss due to final state interactions is predicted to be 
much smaller and exhibit a weaker energy dependence than 
nondissipative hydrodynamics applied throughout the evolu- 
tion. 



There is an ongoing experimental program to produce 
a (transient) deconfined phase of QCD matter Jl 1 in the 
laboratory via nuclear collisions at high energies^ 2] . The 
production mechanism is the liberation of a large num- 
ber of gluons from the nuclear structure functions. The 
plasma is produced from copious minijet gluons at cen- 
tral rapidity, y ~ 0, with transverse momentum px > Po- 
The rapidity density of gluons liberated in central A + A 
collisions can be estimated from [01 



dN 
dy 



(p ) = K T AA (b = 0) J d 2 'p T j ' dx a dx b 



Pt>Po 



s da 



G(x a ,p T )G(x b ,p T ) — -rSfs + t + u) . (1) 
n dt 

s, t, u are the Mandelstam variables for the parton- 
parton scattering process, and da/dt denotes the hard- 
scattering differential cross section in lowest order of per- 
turbative QCD. G(x,p^) denotes the LO gluon distri- 
bution function in the nucleus. The phenomenological 
factor K = 2 accounts approximately for NLO correc- 
tions. The nuclear overlap function Taa(0) = A 2 /-kR\ 
determines the number of binary nucleon collisions in 
head-on reactions within the Glauber approach, where 
Ra — 1.1 A 1 / 3, fm for mass A nuclei. 

For large po, the produced gluon plasma is dilute. As 
p decreases, however, the density of gluons increases 
rapidly due to the increase of G(x,p%) as x « 2pr/y/s 
decreases. It has been conjectured Q that below some 
transverse momentum scale po < p sa t the phase-space 
density of produced gluons may saturate since gg — ► g 
recombination could limit further growth of the struc- 
ture functions. Phenomenologically, this condition may 



arise when gluons (per unit rapidity and transverse area 
7r/p;? at ) become closely packed and fill the available nu- 



clear interaction transverse area. 
p sa t can thus be estimated from 



The saturation scale 



dN 
dy 



(ft 



(2) 



where ~ 1. For 
EKRT H was 



1 the solution reported in 



Psat 



0.208.4 128 
i.34^-0.007. 



■ 0.191 



(3) 



where p sat and y/s are in units of GeV and C\ is the 
average transverse energy per gluon (in units of ftat)- 
The focus of this paper is to investigate whether the 
final observed dEj,/dy can be used to test the pre- 
dicted A and y/s dependence of the initial dE^/dy = 
Cip sa , t dN(p sat )/dy. 

Different gluon saturation models based on non-linear 
evolution and classical Yang-Mills equations (^,^] suggest 
that the factor f3 in (H) may vary parametrically as 



/?(ftat) 



47ra(p sat )A c 
c(A 2 - 1) 



(4) 



where c ~ 1 is a nonperturbative factor proportional to 
the fraction of the initial gluons in the nucleus which 
are liberated. This factor was recently estimated using 
lattice classical Yang-Mills methods || to be c w 1.3. 

The first data || from RHIC on Au+Au collisions at 
y/s = 130A GeV with dN c h/drj w 560 is in fact repro- 
duced by the EKRT saturation model || with 0=1 
assuming isentropic expansion (see On the other 

hand, a fit to the Phobos data using eq.([|) requires 
c w 1.9. Note that the solution of eq. (^) with (3 ^ 1 can 
be obtained from eq. (^|) by rescaling the mass-number 
A — > Aj0 1 l z , and iterating until the stationary point is 
reached. Extrapolating to CERN-LHC energy the mini- 
jet multiplicities are predicted to be dN/dy — 3200 for 
0=1 versus dN/dy — 5100 for /3(p sa t) using eq. (|j). 

While dN(A,y/s)/dy systematics provide one experi- 
mental handle to test different saturation and fixed scale 
models of initial conditions flio| ], another important ob- 
servable that probes collective dynamics is the transverse 
energy per unit rapidity, dEx/dy. In EKRT the final 
value of dE^ /dy was predicted to be much smaller than 
produced initially due to collective longitudinal work as- 
suming the validity of isentropic hydrodynamics. 
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If the expansion proceeds in approximate local equilib- 
rium with pressure p = c 2 e and speed of sound c, then 
the energy density, e(r), must decrease faster than the 
expansion rate T exp — 1/r and leads to a bulk transverse 
energy loss 



E t (t) 
E t (t ) 



r e \tJ 



(5) 



If local equilibrium is maintained during the evolution 
S = cr . In contrast, if the system expands too rapidly 
to maintain local equilibrium, then the effective pressure 
is reduced (relative to that from LQCD) due to dissipa- 
tion. The extreme asymptotically free plasma case corre- 
sponds to free streaming with 6 = 0. Et thus provides an 
important barometric observable that probes the (longi- 
tudinal) pressure in the plasma fill] . There have been of 
course many studies on the magnitude of dissipative ef- 
fects on this and other observables, see e.g. fl2|-|l6l|. The 
new twist on this old problem that we consider here is to 
extend those studies to the novel initial conditions sug- 
gested by gluon saturation models dSlI ■ 

To compute the transverse energy loss due to longitu- 
dinal work, we employ the Boltzmann equation in relax- 
ation time approximation |12|JlqJl7| , 



(6) 



P ■ 9f(p, X) = T re l p ■ U {f eq {p ■ U) - f(p, X)) 



u^ 1 denotes the four-velocity of the comoving frame and 
f eq is the chemical and thermal equilibrium phase space 
distribution, towards which / evolves at a relaxation rate 
T re i. It is important to emphasize that this much sim- 
plified transport equation has been extensively tested 
against full 3+1D covariant parton cascade codes JTJ] 
and provides a surprisingly accurate equation for calcu- 
lating the evolution of the transverse energy observable 
even in highly dissipative systems far from equilibrium 

The relevant relaxation rate is given by the fractional 
energy loss per unit length, 



I-Vei — 



1 dE 
~E~dz~ 



(7) 



which receives a contribution both from elastic and in- 
elastic scattering, 




/167ra 2 A 2 



1 



i, Q 2 

~ l0 S — 2 



(8) 



In these equations p(r) denotes the gluon density in the 
local restframe, p 2 is the Debye screening mass in the 



medium, Q 2 ~ s is the upper bound for the momen- 
tum transfer in the scattering process, and x denotes the 
fraction of energy carried away by radiated gluons. In 
the last step we replaced the momentum transfer q 2 in 
the expression for the radiative energy loss by its aver- 
age, q 2 « p 2 log Q 2 / p 2 . In local thermal equilibrium the 
average energy per gluon and the Debye screening scale 
are 



J/2 = (e/pf * 9T 2 
p 2 = N c g 2 T 2 /3 = AnaT 2 . 



(9) 
(10) 



Assuming that the ratio s/ p 2 is essentially the same even 
out of equilibrium, the relaxation rate is approximately 
given by 



rel 



„ 2 p A L 1 27 
97ra 2 % log- + - T 



P 1 
K in 9na 2 ^log- . 



(11) 



We have set the double- logarithm of Q 2 /p 2 ~ 1/a equal 
to unity. The expression (nTJ) in fact overestimates the re- 
laxation rate at early times, because the screening length 
1//Z. ~ l/gPs&t exceeds formally the horizon at tq = l/p sa ,t 
for longitudinal Bjorken expansion ]^8| and because we 
neglect the suppression of radiation due to formation 
time physics. However, since withp sat < 2 GeV, jw2up 
to the LHC energy domain, we ignore this formal point 
in the discussion below. 

The inelastic, radiative energy loss represents a signif- 
icant source of uncertainty and is especially important 
in chemically undersaturated models of the initial condi- 
tions ]l5| , p^ ]. Within the saturation model, gluon mul- 
tiplication through 2 — > 3 processes may lead to ther- 
malization of the soft radiated gluons at times paramet- 
rically large as compared to To, while the effect on the 
hard part of the gluon distribution is small pp[] . In the 
present paper we do not attempt a more detailed treat- 
ment of radiative energy loss but simply vary the factor 
Kin ~ 1 — 2 to provide a measure of the theoretical un- 
certainites. We note that radiative energy loss for gluons 
with modest pt < 5 GeV, is in any case significantly 
suppressed due to finite kinematic constraints and de- 
structive interference effects ]^T| . 

The EKRT saturation model predicts the gluon density 
to be nearly chemically saturated already at the initial 
time To = 1/psat- This follows from the observation that 
ideal-gas formulas p ~ T 3 , e ~ T 4 applied with chemical 
potential p g — yield the same "temperature" Tq M . 

At the initial time r = l/p aa .t, s/2 = t\j ' p\ = CfpLt- 
Therefore, noting that the comoving gluon density at 
time To is po = p 3 at /ir/3, the ratio of the relaxation rate 
to the expansion rate is given by 



r 



/•< / 



- exp 



= K, 



9a 2 , 1 
~ l0g a 



'/3C 2 



(12) 
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While T re i cx p sat increases as a power of the energy in 
eq. (j^), the Bjorken boundary conditions jl8) force the 
system to expand londitudinally initially also at an in- 
creasing rate T exp (T ) = p sat . The essential quantity 
that fixes the magnitude of the effective pressure rel- 
ative to that predicted by LQCD is the ratio of rates 
in Eq. (|l2]), which dimensionally is simply a function of 
a(Psat)- The asymptotic freedom property |22| of QCD 
therefore requires that this ratio vanishes as y/s — > oo. 
In ([13) the rate of how fast it vanishes is controlled 
by crK in (p sat )/ f3(p sat ). Therefore, with saturation ini- 
tial conditions, asymptotic freedom reduces the effective 
pressure acting at early times r ~ tq and causes the ini- 
tial evolution to deviate from ideal hydrodynamics for 
a time interval that, as we show below, increases with 
energy. 

The total number of interactions during the evolution 
up to time r is given by 



(r) 



Ar'Trei^T 1 ) 



K in 9a 2 log(l/a) 



ft 



log- 



TO 



(13) 



The last expression holds close to the free streaming 
regime ( "Knudsen limit" ) , where the number of scatter- 
ings increases only logarithmically with time. We find 
that <f> reaches on the order of unity at r* 0.4 fm in the 
BNL-RHIC to CERN-LHC energy region. However, the 
local relaxation rate at r* is still less than or on the order 
of the expansion rate (1/r*). In a non- expanding plasma 
(/> ~ 1 provides a rough equilibration criterion. However, 
as long as r re ; is not significantly larger than T exp this 
criterion is insufficient to address how much collective 
hydrodynamic work can be done by the plasma. 

For a quantitative estimate, we must solve the kinetic 
equations (^J). The first energy moment of that equa- 
tion together with energy conservation to replace e eq by 
e(r), results in the energy density evolution equation 
0,f6|JT7| 



TC 



f >') 



"(0) 



(14) 



This equation applies if the initially produced partons 
have a vanishing longitudinal momentum spread in the 
comoving frame, i.e. assuming a strong correlation be- 
tween space-time rapidity and momentum space rapid- 
ity p8[ . The function h(x) appearing in jl^ ) is given by 
2h(x)VT 



x" = x\fl — x'- + arcsm V 1 — x* 
that S — > 1/3 as t - > oo . 

Expanding ( |l4| ) to first order in (f> yields 



it insures 



re 
Toeo 



giving for r ~ tq 



7T 

To 



\ 9K, 



ftC\ 
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log- 
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(15) 



(16) 



Note that S — > as p sa ± — > oo in accordance with the 
discussion above. 

For our numerical estimates we let a creep with time. 
Since the effective temperature scales as e(r) 1 / 4 , the ef- 
fective coupling, a(T e ff), increases slowly with time ap- 
proximately as 



a(r) 



QCD 



1/2N 



, (17) 



with Aqcd — 200 MeV. Furthermore, longitudinal ex- 
pansion constrains p(t)t to remain constant. As the sys- 
tem cools, the mean center of mass energy in collisions 
also decreases as 



s(r) 



P(r) 



s(r ) 



Te(r) 
to£(to) 



(18) 



The number of collisions between tq and r in this case is 
given by 



r, e i(r ) / dr' . 

Jro V a ( T 0) 



a(r')\ 2 fe(r o y 2 



e(r') 



X 



log(l + l/ajT 1 )) 
log(l + l/a(r )) 



(19) 



We regulated the logarithmic dependence above for nu- 
merical stability. Note that with eq. (pTj|), eq. (14) is a 
nonlinear self-consistency equation for e(r). 
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FIG. 1. The ratio of effective longitudinal pressure to 
energy density as a function of the energy density along the 
dynamical path is shown for SPS, RHIC, and LHC saturation 
initial conditions 

|). Solid (dashed) curves are for (3 = 1 and 
K in = 1(2). The LQCD equation of state |l) is also shown 
for comparison. 

In Figure [l] we show the effective longitudinal pressure 
as a function of the energy density for yfs = 20, 200, 5400 



3 



A GeV saturation initial conditions. The ratio p/e is 
defined here by S(t) as obtained solving eq. ( |T4| ) numer- 
ically. Our definition of the effective pressure absorbes 
all disspative corrections to the perfect fluid equation, 
7V<9„[(e + p)u^u" - pg^ v \ = de/dr + (e + p)/r = 0. 
For comparison the pressure of equilibrium QCD is also 
shown for Nf = 3. This curve is obtained from the 
Nj = lattice data of fl rescaling the number of rel- 
ativistic degrees of freedom by 47.5/16, and assuming a 
transition temperature T c (Nf = 3) = 160 MeV. 

Initially p/e starts at zero in this model and remains 
small for a large time relative to l/p s &t because the 
plasma is torn apart by the initial rapid longitudinal ex- 
pansion. The effective pressure approaches the LQCD 
curve from below and reaches it at a time tl ~ 1 — 2 fin 
at RHIC y/s = 200 AGeV, by which time the energy den- 
sity has dropped by an order of magnitude, = e(rj,) = 
6.5 - 12 GeV/fm 3 . For LHC y/s = 5400 AGeV, t l = 
3 — 7 fm during which the energy density falls by almost 
two orders of magnitude to e£ = 9.5 — 21.5 GeV/fm 3 . 
The quoted intervals correspond to Ki n = 1 — 2 using 
the EKRT parametrization (Q). The contrast between 
the dynamical path followed by the saturated plasma 
compared to the equilibrium equation of state is striking. 
A qualitatively similar behavior of the early logitudinal 
pressure has also been found from solutions of diffusion 
equations E3). 




0.1 1/2 1.0 
s /A (TeV) 



10.0 



FIG. 2. The ratios of the final to the initial transverse en- 
ergy per unit rapidity are shown as a function of beam energy 
for central Au+Au collisions. The initial value corresponds to 
the EKRT parametrization (f3 — 1) and the transport results 
are for Ki„ — 1,2. For comparison, the final transverse energy 
assuming that local equilibrium was maintained throughout 
the evolution is also shown. 

Since the longitudinal gradients at tl are much smaller 
than at To, the evolution beyond tl is much more likely 
to follow isentropic hydrodynamics along the lattice QCD 
equation of state. In this case one could calculate detailed 



differential hadronic observables along the same lines as 
in |24j ] using the conditions at Tl as the initial conditions 
for 3+ ID hydrodynamics. 

The main experimentally observable consequence of 
the reduced effective pressure is shown in Fig. 0. The ra- 
tio Et/Et — Tf€f /roeo has been obtained from the solu- 
tion of the transport equation assuming e/ = 2 GeV/fm 3 
which corresponds roughly to T ~ T c . Tf is estimated 
assuming hydrodynamic expansion from the point were 
the trajectories in Fig. |l| reach the LQCD curve at time 
tl , with the equation of state p/e = a + b log e (e in units 
of GeV/fm 3 ). The parameters a = 0.051, b = 0.092 pro- 
vide a reasonable fit to the LQCD curve shown in Fig. [I]. 
In this case, 



TL 

TL 



b log e L 



1 



b log e / 



1/6 



(20) 



On the other hand, if ideal hydrodynamics were appli- 
cable already at To, the final observed transverse energy 
for 1+1 dimensional adiabatic expansion would be 



f 



T f£f = T f( T f s f-Pf) _ T± 
T e T (T s —po) T 



(21) 



The last step follows both for poj ~ as well as 
Poj = TqjSqj /4 from the condition of entropy conser- 
vation, ts = const H. Strong transverse expansion leads 
to slightly larger Elp but we shall neglect that small effect 
here for simplicity. Clearly, for Tf ~ T c ss 160 MeV one 
would observe a much smaller transverse energy in the 
final state than in the initial state. Moreover, EL/Ej, 
would also have significantly stronger energy dependence 
such that E T deviates more and more from E l T with 
increasing y/s. In this sense isentropic hydrodynamics 
erases information on the interesting initial conditions via 
this observable. The solutions of the transport equations 
clearly show a smaller decrease of Ej, and of the logarith- 
mic slope, k — dlogi?^/dlog\/s, due to final state in- 
teractions. We find that k = 0.50 for the initial state (||) 
evolved with Ki n = 1, k = 0.46 with Ki n = 2, while 
K = 0.40 with isentropic expansion, eq. (pl|). For compar- 
ison, the initial EKRT saturated E l T — irR^Toeo scales 
with the higher power k = 0.59 according to eq. (|3|). The 
fractional transverse energy loss is thus less dependent on 
energy than for entropy conserving expansion for which 
Ej\/Et cx 1/Xo oc 1/^/s ' 2 . This is due to the increas- 
ingly long time spent far from equilibrium in Fig. |l| as 
the beam energy increases. 

The results in Fig. |^ are encouraging from the point 
of view of searching for evidence of gluon saturation in 
nuclei at high energies. Experimental data on dEx/dy 
or dEx/drj for central Au+Au collisions at RHIC will 
soon provide a new test of saturation and non-saturation 
models at those energies. Since we predict that dissipa- 
tive effects reduce considerably the effective longitudinal 
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pressure in Fig. 1, the beam energy dependence of the 
transverse energy is expected to reflect much more accu- 
rately the predicted power law dependence of the initial 
conditions as seen in Fig. 2. We therefore conclude that 
the energy and A systematics of the bulk calorimetric 
observable, dE'r/dy, will be a sensitive test of saturation 
models of gluon plasmas produced in the RHIC to LHC 
energy range. 
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